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For a model of DNA denaturation, exponents describing the distributions of denaturated loops 
and unzipped end-segments are determined by exact enumeration and by Monte Carlo simulations in 
two and three dimensions. The loop distributions are consistent with first order thermal denaturation 
in both cases. Results for end-segments show a coexistence of two distinct power laws in the relative 
distributions, which is not foreseen by a recent approach in which DNA is treated as a homogeneous 
network of linear polymer segments. This unexpected feature, and the discrepancies with such an 
approach, are explained in terms of a refined scaling picture in which a precise distinction is made 
between network branches representing single stranded and effective double stranded segments. 
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I. INTRODUCTION 



A DNA molecule may undergo transitions from a dou- 
ble stranded to a single stranded state, either under the 
effect of an increase in temperature T (thermal denatu- 
ration) , or through applied forces at one end of the chain 
(mechanical unzipping) In the characterization of 

such transitions, and in the determination of their uni- 
versal, asymptotic features, substantial progresses were 
made recently by applications of models and methods of 
polymer statistics 1^-^. Among these progresses is an 
extension |^ of the classical Poland and Sheraga (PS) 
model In the PS model the partition function of a 
DNA chain is approximated by that of a sequence of non- 
interacting double-stranded segments and denaturated 
loops, and the thermal denaturation transition results 
of second order type Recently, excluded volume 

effects between a loop and the rest of the chain were 
included in the PS description in an approximate way 
[pI, using results from the theory of polymer networks 
[p|. This approach predicts a first order denaturation, in 
agreement with very recent numerical studies of models 
taking fully into account the self and mutual avoidance 
among loops and double segments 0,11 . Quite remark- 
ably, the approximate scheme of Ref. ||^ yields results 
which are in good quantitative agreement with Monte 
Carlo simulations Q. Most recently, predictions based 
on the theory of polymer networks were also made for 
the case of mechanical unzipping . 

Besides confirming the expected first order character 
of thermal denaturation, the results of Ref. |^ demon- 
strated that excluded volume effects alone are respon- 
sible for this character, while the difference in stiffness 
between double and single stranded DNA, and sequence 
heterogeneity do not affect the asymptotic nature of the 
transition. 



All these results rise interesting and debated |11| is- 
sues and open new perspectives in the field. First of all, 
one would like to test numerically the existing analyti- 
cal estimates, in particular the new ones pertaining to 
the case of mechanical unzipping. The validity of crucial 
predictions relies on such tests, which should also reveal 
whether the quantitative success in the case of thermal 
unzipping is rather fortuitous, or there is some deep and 
systematic basis for it. The network picture proposed 
in Ref. besides allowing some elegant and successful 
estimates, could constitute an important step forward in 
our way of representing the physics of DNA at denat- 
uration. Progress in the assessment of the validity and 
limitations of a network picture for denaturating DNA 
should be allowed by a more careful and systematic anal- 
ysis of numerical results for specific models. The possible 
extension and improvement of previous analyses for the 
relevant three dimensional case is also a main motivation 
of the present work. 

In this article we consider a lattice model of DNA both 
in two and in three dimensions (d) with excluded volume 
effects fully implemented. By various numerical methods 
we estimate length distributions for denaturated loops 
and unzipped end-segments. In particular, for the transi- 
tion in three dimensions we give here exponent estimates 
which extend and improve the results of a previous study 

. Moreover, our results allow to address the basic issues 
mentioned above, and to clearly identify some qualita- 
tive and quantitative limitations of the picture proposed 
in Ref. . The analysis gives also hints which allow 
us to propose a generalization of the polymer network 
representation of denaturating DNA. Our model in two 
dimensions, besides deserving some interest in connec- 
tion with problems like the unzipping of double stranded 
polymers adsorbed on a substrate, offers an ideal context 
in which to compare the predictions of Refs. PpOl with 
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numerical results. Indeed, for d = 2 those predictions 
are based on exactly known network exponents, while 
in d = 3 the same exponents have been approximately 
determined. 

The model studied here was introduced in Ref. and 
further analyzed and extended in Ref. ||^. We consider 
two self-avoiding walks (SAWs) of length N on square 
and cubic lattices, described by the vectors identifying 
the positions of each monomer fi{i) and ^2(1), with 
< i < N. The walks, which represent the two DNA 
strands, have a common origin ri(0) = f2(0) and only 
monomers with the same coordinate along the SAWs can 
overlap each other (i.e. ri{j) = r2{k) only if j = k). An 
overlap corresponds to a bound state of complementary 
DNA base pairs, to which we assign an energy e — —1, 
thus neglecting the effects of sequence heterogeneities. 
This is an acceptable approximation in view of the fact 
that, at coarse grained level, each monomer (site vis- 
ited by the SAW) should represent a whole persistence 
length of the single strand, which includes several (« 10) 
bases. A denaturated loop of length / occurs whenever, 
for some i, r2{i) = ri{i), r2(* + ' + 1) = + I + 1), and 
r2{i + k) ^fi{i + k), fc = 1,2,..,L 

At a temperature T each configuration lu of the two 
strands appears in the statistics with probability propor- 
tional to its Boltzmann's weight exp{— (3 H (uj)} , where 
f3 = 1/T (Boltzmann's constant = 1) and H is e times 
the number of bound base pairs in uj. The resulting be- 
havior of the DNA model resembles the scaling of a 2A^- 
step SAW as long as the inverse temperature f3 is lower 
than a critical value /3c. For (3 > Pc the strands are typi- 
cally paired in a sequence of bound segments alternating 
with denatured loops, the latter gradually shrinking and 
becoming rarer for increasing /3. For f3 > (3c and iV — » 00 
the scaling of an iV-step SAW holds. Around f3c there is 
a crossover from one regime to the other, and exactly at 
f3c peculiar scaling behaviors are expected for the distri- 
butions of loops and end-segments. 

II. THE TWO DIMENSIONAL CASE 

We focus first on the two-dimensional case. Fig. plots 
the logarithm of P{1, N), the pdf of finding a denaturated 
loop of length / within a chain of total length iV, as a func- 
tion of In/ at the transition point. The behavior of this 
pdf at denaturation determines the order of the transi- 
tion The data are obtained from exact enumera- 
tions for walks up to iV = 15. For each chain length, the 
pdf is sampled at the temperature corresponding to the 
specific heat maximum. This choice offers the advantage 
that, if the specific heat diverges at the transition, like 
in our case, the correct transition temperature is asymp- 
totically singled out for N 00. Thus, in this limit, the 
sequence of distributions should automatically approach 
the critical pdf, for which we expect power-law scaling. 




FIG. 1. In P{l,N) vs. Inl for various N values. Upper in- 
set: In P{N/2, N) vs. InA^; data are well-fitted by a line from 
which we estimate c = 2.44(6). Lower inset: Scaling collapse 
of the data. 




FIG. 2. (a) EfTective exponent 1/ of the end-to-end dis- 
tance, for e A = {80, 120, 160, 240} (dashed curve), 
N € B = {160, 240, 320, 480} (dotted curve), and iV e AuB 
(continuous curve), as a function of /3 in d = 2. The horizon- 
tal dotted line marks the exactly known SAW = 3/4, while 
the vertical line is our estimate f3c = 0.7525(3) (the error 
is indicated by the gray band), (b) Similar plots for d — 3 
{A = {80, 120, 160, 240}, B = {120, 160, 240, 320}). In this 
case the intersections are even better localized around the 
expected SAW v « 0.5877. The transition temperature de- 
termined in this way is almost coincident with the estimate 
/?c = 1.3413(4) |] indicated by the vertical line and the gray 
band. 



Finite-size effects should be described by P{1, N) = 
l'^q{l/N), with q a suitable scaling function. In order to 
determine c, we consider, e.g., P{N/2, N), which should 
scale cx N^'^ for N 00. Such quantity is shown in 
the upper inset of Fig. |l|, plotted as a function of N in 
a log-log scale. From a linear fit of the data we obtain 
the estimate c = 2.44(6). The lower inset shows a scal- 
ing collapse of the data, obtained with c = 2.44. The 
good quality of the collapse indicates that the system is 
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very close to the asymptotic regime although the chains 
are quite short. We also performed a Monte Carlo de- 
termination of c using the pruned enriched Rosenbluth 
Method (PERM) through which walks are gener- 

ated by a growth procedure. In this case the critical 
temperature at which the pdf was sampled was deter- 
mined by carefully monitoring the scaling with N of the 
average end-to-end distance {\r2{N) — 'ri{N)\) as a func- 
tion of temperature. The effective exponents describing 
the growth with TV of this quantity are plotted as a func- 
tion of (3 for different chain lengths in Fig. || The 
intersections of the various curves in a very narrow range 
signals the crossover expected at denaturation for this 
quantity and allows to locate the melting temperature 
rather accurately, i.e. (3c = 0.7525(3). By fitting the ini- 
tial slope of the critical pdf for chains up to iV = 480 
we extrapolate c = 2.46(9), in good agreement with the 
exact enumeration results. 




FIG. 3. (a) Configuration of a polymer network with £, — 2, 
ni — 5, ns = 1 and 714 = 3. (b) Loop of length / embedded in 
a chain of length N — I and (c) chain of length N — m with 
bifurcating end segments of length m. 

The value of c for loops with I <C N was predicted an- 
alytically ||] using exact results for cntropic exponents of 
networks of arbitrary topology For a network of 

fixed topology G (see example in Fig. ^(a)), with n self- 
and mutually avoiding segments, Duplantier on the 
basis of renormalization group arguments, postulated the 
following scaling form for the total number of configura- 
tions: 

TG-^i^N'^^-' ff^l/-!,../-^) (1) 

where k is the length of the z-th segment, N — k and 
/ a scaling function. The value of 70 depends on the 
number of independent loops, £, and on the number of 
vertices with k legs, rifc, as 

7g = 1 - Cdv + ^ rifcCTfe (2) 

k 

where v is the radius of gyration exponent and ak, 
k = 1,2, .., are exactly known exponents in d = 2 M: 



(2 - k){9k + 2) 



This general scaling framework was applied to the 
DNA unzipping by considering relevant network topolo- 
gies for the problem For example, in order to study 
the denaturated loop length pdf, one can assimilate the 
situation of a typical loop within DNA to that of the loop 
in Fig. |^(b). This amounts to assume that the action of 
the rest of the DNA molecule is the same as that of two 
long linear tails, thus totally disregarding the presence of 
other loops. According to Eq. (|l|), for G corresponding 
to the topology in Fig. ||(b), one has j|] 

Tioop ~ ^i^^N-'^•'+^'''+^''^h{l/N) . (4) 

In this equation we explicitly assume that, while the con- 
nectivity of a loop step is equal to /i, that of a step of the 
tails is /i^. This assumption is consistent with the ther- 
modynamics of the denaturation transition and follows 
from the continuity of the canonical free energy of the 
system, and from the fact that each step of the tails cor- 
responds in fact to two steps of the loops. In this way the 
I dependence of the r.h.s. of Eq. (||) does not enter in the 
exponential growth factor, an important requisite for the 
derivations below |l5j . For N ^ I the number of configu- 
rations should reduce to that of a single double stranded 
chain of length N, i.e. F - fj.'^^N'^-'^ with 7 = 1 + 2cri, 
according to Eq. (|l|). This requires h{x) x^'^^^'^"^ , for 
X 1. So, the / dependence in Eq. (^) becomes ~ 
with (I 

c = dv — 2(73 ■ (5) 

Clearly, in this approximation, c is also the exponent by 
which the loop length pdf, P, scales. In reality, the seg- 
ments departing from the two sides of the loop replace 
more complex fluctuating structures containing denatu- 
rated bubbles of all sizes, separated by short linear double 
stranded segments. In d = 2, v = 3/4 and 173 = —29/64, 
therefore c — 2 + 13/32 w 2.41, which is a value consis- 
tent within error bars with our numerical estimates. As 
already known in d = 3 ||], this agreement implies that 
the sequences of loops which "dress" the two segments 
departing from the loop in Fig. ^(b) have very little ef- 
fects on the value of c. 

We consider now the distribution of end-segments. 
With the assumed boundary conditions, denaturated 
end segments of length m occur in configurations where 
'ri{N — m) — f2{N — m) while ri(fc) ^ r2{k) for k > 
N — m. The statistical geometry of denaturated end- 
segments is expected to be relevant for situations occur- 
ring in mechanical unzipping experiments pO| ]. Indeed, 
as a rule, this unzipping is induced by applying forces 
which separate the strand extremes fi{N) and f2{N) by 
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micromanipulation techniques p6| . In the same spirit as 
in the case of denaturated loops, one can consider now 
the network geometry of Fig. 0(c) un , for which 



r 



fork 



(6) 



By matching again, for m <C iV, with the partition func- 
tion of a SAW of length N, one gets for the distribution 
of the end segment lengths Peim, N) ^ m~'^ with 



- (O-I + (73 ) 



(7) 



One should note that this last result holds also in case 
one tries to match the behavior of Ffork with that of a 
simple SAW of length 2N in the hmit (iV - m) < N. In 
this case the pdf results ^ {N — m)~^. 




FIG. 4. Plot of In Peim, N) vs. Inm (a) and ln(7V - m)(b) 
at the estimated critical point /3c = 0.7525 and for — 80, 
120, 240, 320 and 480 (from top to bottom). 

Figure || shows a plot of the logarithm of the pdf 
for end segments of length m, lnPe(TO, A^), as a func- 
tion of (a) Inm and (b) ln(A^ — to). The data ob- 
tained from the PERM show coexistence of two dis- 
tinct power-law scaling behaviors, i.e. Pe{m,N) ~ m~''^ 
and Pe{m, N) ^ (A^ — m)^'^^ where a linear fit yields 
ci = 0.23(1) and C2 — 0.35(1), respectively. The slope 
determinations are rather sharp in this case, and, in fact, 
the criterion of selecting the denaturation temperature 
in correspondence with the simultaneous manifestation 
of the power law behaviors of Pe is very efficient, and 
consistent with that based on monitoring the behavior of 

{\UN)-n{N)\). 

The existence of two distinct slopes, and the values 
of the exponents are in disagreement with what one ex- 
pects on the basis of the network approximation, i.e. 
c = 9/32 w 0.28 in d = 2. Thus, for some reason, 
on Pe{m,N) the effects of the structure of the double 
stranded part of the chain are noticeable and the schema- 
tization through a simple polymer network topology is 



not fully adequate to represent the physics. In the case 
N — m <^ N the number of single chain configurations is 
indeed that of a simple linear SAW of length 2A^, which 
is asymptotically exactly known in two dimensions [Q. 
On the other hand, for to <C A^, the configurations to 
count are those of an effective linear SAW chain of length 
A^, whose internal structure contains loops at all scales. 
It is conceivable that not only the connective constants, 
but also the entropic scaling properties of such effective 
walk differ from those of a standard SAW. The difference 
must concern the power law correction to the exponential 
growth factor of the number of configurations as a func- 
tion of A^. Indeed, as already remarked above, the con- 
nectivity constant of the effective, double stranded walk 
at the transition, must be just the square of the simple 
SAW connectivity constant. This follows from an obvious 
requirement of free energy continuity and from the fact 
that in the high temperature region the two strands be- 
have as unbound simple SAW, with a total length twice 
that of the effective double stranded walk, independent 
of temperature ||l^ . 

To investigate this issue further we determined directly, 
on the basis of PERM data, also the overall entropic 
behavior of the DNA chain. We considered two types 
of boundary conditions: (1) the extremes (fiiN) and 
■r2{N)) of the chain are free and (2) are forced to join in 
a single point {ri{N) — r2(A^)), while in both cases, the 
strands have still a common origin (ri(0) = '^2(0)). 

Condition (1) is the one applying to the effective walk 
discussed above. We indicate with Z^-* and Z^' the 
corresponding partition functions. In the spirit of the 
network approximation, we neglect the contribution of 
denaturated loops within the double stranded phase one 



has Z 



(2) 
N 



2N 



.(1) 



A^T \ with 7(2) = 1 + 2(71. We esti- 
mate also Z^j^' within the same general framework, by 
integrating the partition function Ffork of Eq. (^) over all 
possible end-segment lengths. This integration gives 



dm F 



fork 



2Ar^3cri+cr3+l 



(8) 



as the integration of the scaling function g{m/N) yields 
an extra factor proportional to N . Defining zl!^ ~ 



2N 



one eventually gets: 
7(1) = 2 + 3(71 4 



f73 . (9) 

This is a theoretical expression for the entropic expo- 
nent of a DNA molecule as a whole, and can be directly 
compared with numerical estimates, one of which already 
exists in d = 3 Ifzl , as we discuss in the next sections. 

/2) 



To calculate the entropic exponents 7^^^ and 7^ 



we 



estimated the quantity (Z^'^f/^^V^^^ by PERM sam- 
pling for reasonably long chains. For N 00 one expects 



7(0 ' 

4\ 



1 + In 2 



2N 



i = 1, 2. (10) 
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Since /z must coincide with the SAW connective con- 
stant, which in d = 2 is very precisely known (/x = 
2.63815852927(1) ||l|]), large TV data for the quantity on 
the l.h.s. of Eq. ( |10[ ) can be fitted by keeping as unique 
fitting parameter 7^'^ in the r.h.s. expression. Alterna- 
tively, one can assume a value for 7^*^ and check whether 
data appear consistent with the assumed correction term 
in the same r.h.s. expression. This consistency test is 
best applied to sets of data pertaining to different tem- 
peratures close to the transition, as illustrated in the next 
section. Since the correction term is rather sensitive to 
the choice of the temperature, this offers also a way of 
locating the transition. Here we just mention that the 
slopes predicted on the basis of the network approxima- 
tion in d = 2 appear slightly, but definitely, inconsistent 
with the data, as discussed in the next section. 

III. DENATURATING DNA AS A COPOLYMER 
NETWORK 

The polymer network representation of denaturat- 
ing DNA was originally proposed as a useful, but 
definitely approximate tool, without the possibility of 
controllable and systematic improvements. For exam- 
ple, the environment seen by a single fluctuating loop 
within the molecule was not proposed as something 
unique, and choices slightly different from that discussed 
above were also discussed ||^,|l^. With these alternative 
choices, leading to slightly different results, the portions 
of the molecule surrounding the loop were not necessarily 
treated as made of simple double segments. 

Even if the approximate character remains, all the dis- 
crepancies and inconsistencies discussed in the previous 
section can be resolved by a refinement of the whole pic- 
ture and an improvement of the approach, which could 
allow more accurate predictions in the future. In the 
new perspective we propose, the rules for associating a 
network schematization to loops or end-segments should 
be regarded as unique. This follows from the fact that 
the new polymeric entities one defines are supposed to 
account, at an effective level, for all the complications 
arising from the fluctuating geometry of the model. 

The legs of the network are either simple self- avoiding 
chains, in case they are really made of single strands, 
or effective, dressed segments, when they represent fluc- 
tuating double stranded portions of the molecule com- 
prised between two given bound base pairs. This clear 
cut distinction suggests that the exponents ct^, associ- 
ated with the network vertices, should be modifled with 
respect to the "bare" values considered so far, as soon 
as at least one outgoing leg is of the double stranded, 
dressed type. Indeed, within the framework of a contin- 
uum Edwards model description of the inhomogeneous 
polymer network, and of renormalization group ideas, 
one can postulate the validity of homogeneity laws and 



exponent relations analogous to Eqs. (|^) and (||), with 
modified cr's where appropriate. Something like this was 
already done in Ref. ||l^ where explicit field theoretical 
calculations were performed in the case of copolymer star 
networks (no loops). These stars were composed of mu- 
tually excluding branches made either of random walks, 
or of self avoiding walks. In the case considered here the 
copolymer is made of two kinds of segments which do not 
differ as far as metric scaling exponents are concerned 
(same v) , but can have different entropic scaling proper- 
ties. It is indeed quite conceivable that an effective linear 
structure, which should be resolved into a sequence of de- 
naturated segments and loops of all sizes, could have an 
entropic scaling different from the one of a simple SAW 
on the lattice, even if the elongation grows in the same 
way with N. 




FIG. 5. Representation of the geometries relevant for the 
DNA problem by copolymer networks. Thick lines denote 
branches of the network composed by an alternating sequence 
of bound segments and denaturated loops, while thin lines are 
genuine single stranded segments. Only the latter have the 
same metric and entropic scaling properties as a SAW. 

The above considerations lead us to introduce two new 
exponents: ai = ai + Aai for an isolated vertex with a 
dressed segment, and W3 — + Aa^ for a vertex join- 
ing a dressed segment and two self avoiding walks, as 
schematically illustrated in Fig. ||. As discussed above, 
the introduction of such ct's amounts to postulating a 
generalization of the scaling in Eq. (|^) to copolymer net- 
works with two different types of segments. Following 
the same arguments leading to Eqs.(||), (^ and (||) and 
using now the modified exponents wherever appropriate, 
one finds: 



Cl = - (o-i + (T3) - Acts + Acti , (11) 

C2 = - {(71 + as) - Acts - Aai , (12) 

c = di^ - 2(T3 - 2Acr3 , (13) 

7(2) = 1 + 2(71 + 2A(Ti , (14) 

7^^^ =2 + 3ai+a3 + Acti + Acrg , (15) 



which is a set of five equations with only two unknown 
parameters, Aai and Aa2- These equations should be 
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regarded as consistency requirements in order to test 
the validity of the proposed copolymer picture. Notice 
that within the copolymer network scheme the expo- 
nents associated to the distributions of short and long 
end-segments are distinct, ci ^ C2, as soon as Acti ^ 0. 
By solving the first two equations (|l^) and with the 
numerical values for c\ and C2, we find Atra — —0.01(1) 
and Ai7i = —0.06(1). While the former is actually com- 
patible with zero, the latter is not. Once the values of 
Aci and Aca have been fixed, we can check for the con- 
sistency of the other exponents using Eqs. (p^, (14) 
and (p^). Inserting the calculated Acts into Eq. (12) 
we find c = 2.43(2), i.e. a value slightly higher than 
what predicted from Eq. (^). Indeed, our numerical esti- 
mates suggest for c a slightly higher value than that pre- 
dicted on the basis of the bare network approximation, 
although the error bars cover both values. Next we con- 
sider the last two equations (^J) and (|l5|); as both A(7i 
and A(T3 are negative we expect that 7'^and 7^^^ would 
be somewhat smaller than the values predicted from the 
homogeneous network approximation. Substituting the 
above numerical values of Acii and Acra into Eqs. (|lj) 
and (|l|) we find 7(2) « 1.22(2) and 7(1) « 1.99(2) 
(recall that the homogeneous network predictions are 
7(2) = 1 -Ml/32 « 1.34 and7(^) = 2 1/16 « 2.06). We 
first examine case (1). Since the differences between the 
"bare" and "dressed" values of 7^^-' are small, to magnify 
the asymptotic details we considered the quantity 



We conclude that the improved representation based 
on the idea of a copolymer network with modified en- 
tropic exponents allows to match well all our numerical 
results, with just two adjustable parameters. 
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FIG. 6. The function /(A^), defined in Eq. (|l6|), for 
N 00 should approach fi with zero slope if the 7^^' ex- 
ponent would be correctly determined by the homogeneous 
network approximation. Instead, the trend predicted by the 
copolymer theory (continuous line) fits rather well the data 
(/3 = 0.7520, 0.7523, 0.7526, 0.7529, 0.7532, from below. We 
recall that we estimate /3c = 0.7525(3)). 



IV. THE THREE DIMENSIONAL CASE 



with 7* = 33/16. The coefficient of the 1/N term in the 
r.h.s. of the previous expression has been chosen such 
that if z'j^^ would scale with the homogeneous network 
exponent then f{N) should approach the connectivity 
constant fi with zero slope when plotted as function of 
1/2N (as follows from Eq. ([l0|)). A 7'-'^' larger (smaller) 
than its homogeneous network value would imply a f{N) 
approaching fi with a positive (negative) slope, equal to 
(^(i)_^*)^ln2. Figure I shows a plot of /(TV) vs. 1/2N 
for five different temperatures around the estimated criti- 
cal one. Clearly f{N) approaches /i in the limit of large N 
with a negative slope. The solid line represents the slope 
for the value of 7*^^^ = 1.99 as obtained from Eq. (15), 

which apparently fits reasonably well the data. 

(2) 

In the case of Z^' , from all configurations generated 
by the PERM only those where the two strands have 
common endpoints {fi{N) — r2iN)) are considered and 
this reduces considerably the statistics available. The nu- 
merical results we find are consistent with both proposed 
pictures, since the low precision of data, due to the in- 
sufficient sampling, produces large statistical errors. So, 
the case (2), while also giving meaningful data, does not 
help in the determination of the right scenario. 



We now consider the d = 3 case, where homogeneous 
network exponents can be deduced from estimates of 
7=1 + 2(Ti |20| (cti 0.079) or from e-expansion results 
combined with resummation techniques (cts ~ —0.175) 
|l|- From Eqs.d) and @ one gets c « 2.11 and 
c « 0.095. Previous Monte Carlo simulations yielded 
c = 2.10(4) In the present work we made an ex- 
tra effort in order to get a reliable estimate of the var- 
ious exponents. With an analysis of the scaling behav- 
ior of (|ri(A^) — f2{N)\) we get here a very precise es- 
timate of the melting temperature [see Fig. ||(b)], con- 
sistent with that of Ref. . Extensive PERM sampling 
of the loop distribution at the estimated transition tem- 
perature yields c = 2.18(6) (see Fig. |^). This refined 
value, confirming the first order character of the tran- 
sition, is slightly higher than that reported in Ref. 
although compatible within the uncertainties. The dis- 
crepancy could be imputed to a slight overestimation of 
the transition temperature made in that reference. An 
analysis of the end-segments distribution yields, as in 
d = 2, two slightly different exponents, C2 > ci, signalling 
again deviation from the prediction of the homogeneous 
network approximation. We estimate C2 = 0.16(1) and 
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ci = 0.14(1). These estimates are rather sharp, and im- 
ply AcTi = 0.01(1) and Acrg = 0.055(10). By insert- 
ing these values for the Act's into our expression for c, 
Eq. (pi]), we get c « 2.22, which is well compatible with 
the last mentioned estimate illustrated in Fig. M. 




1 2 3 4 5 6 

ln(0 



FIG. 7. Log-log plot of the loops pdf at the critical point 
as function of their length for chains of various lengths. 



d = 2 case, we set 7* equal to the value derived within 
the homogeneous network picture for 7«, 7* = 2.06). 
The value of the connectivity constant is ^ = 4.68404(9) 
for the cubic lattice a value indeed approached by 
all data sets for (3 < (3c- At high temperatures the 
data clearly show SAW scaling, as expected, while at 
the transition point f{N) seems indeed to approach a 
slope given by the "dressed" exponent 7^^^ — 2.00 (solid 
line). To confirm this, at the estimated critical point, we 
performed a series of calculations up to very long chains 
{N = 2000), but with somewhat lower statistics. The 
latter data are plotted in gray in Fig. ^: indeed they 
seem to follow quite closely the solid line, as predicted 
by the inhomogeneous network theory. 

Thus, even if the error bars on some exponents are 
still relatively large, the corrections suggested by a block 
copolymer network representation of DNA seem to be 
well consistent with the numerical scenario, like in the 
two-dimensional case. 
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FIG. 8. As in Fig. |: from below, (3 = 0.293 (SAW 
regime, evidenced by the dotted line, i.e. 7* — 7, with 
gamma representing the SAW exponent in ci = 3 pc|]), 
1.335, 1.3407, 1.3413(/3c), 1.3419. The gray band shows data 
at /3c and longer A'^ (up to 2000), but with lower statistics. 

Causo et al. determined in d = 3 the exponent 
7^^) ~ 2.09(10), a value which should be compared with 
the Eqs. (|) and (^ derived in this paper. Inserting the 
above determined values for Aai and Acts in Eq. (|l5| ) 
we obtain 7^^-' « 2.00, to be compared with a value 
7^^) K, 2.06 if one sticks to the homogeneous network 
prediction, given by Eq. (^. As the 7*^^^ of Ref. is 
compatible with both above values we made an effort to 
estimate it again with improved accuracy. 

Figure |^ shows a plot of f{N) vs. 1/2N for five dif- 
ferent temperatures of which three around /3c (as for the 



V. CONCLUSION 

We investigated a lattice model of DNA denaturation 
both in d = 2 and d = 3 and determined the exponents 
associated to the decays of the pdf 's for loops and end- 
segments. For the loop pdf exponent we find c > 2, 
which implies a first order denaturation as the average 
loop length remains finite at the transition point . In 
d = 2 the first order character is even more pronounced 
than in c? = 3 and the values of c are quite consistent 
with those predicted analytically on the basis of entropic 
exponent of homogeneous networks An analysis of 
the end-segment lengths reveals that the corresponding 
pdf displays two distinct power law behaviors, one apply- 
ing at small, and the other at large fork openings. This 
unexpected feature, not predicted within the framework 
of the approximate description based on homogeneous 
polymer networks |6| , suggests to describe the DNA fluc- 
tuating geometry as a copolymer network, in which a 
distinction is made between single stranded and effec- 
tive double stranded segments. The latter are assumed 
to have different entropic exponents than SAW's, which 
lead to the introduction of a generalization of the previ- 
ously known expression for pdf and entropic exponents 
(see Eqs. ( p[ - |l5| )). We tested the compatibility of the ob- 
served results with this effective copolymer network pic- 
ture, which is assumed to catch the essential physics of 
the transition. Within this new framework, the different 
scaling behaviors of the end-segment pdf can be qualita- 
tively explained, and the notion of overall entropic scaling 
of the macromolecule acquires a precise and consistent 
meaning. Especially in d ~ 2, the numerical evidence 
that a block copolymer network picture is compatible 
with the overall data, is rather clear and convincing. 
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A main reason why the copolymer network descrip- 
tion turns out to be well compatible with the observed 
scalings is probably the fact that the pdf of loop length 
scales always with a sufficiently large c exponent. So, 
even for N — > oo, the average width of the loops, and 
thus also of the dressed segments, remains finite. This 
could be an important requisite for the validity of the 
network picture. Indeed, it would be interesting to test 
whether similar copolymer pictures work also for other 
unzipping transitions of double stranded polymers with 
a smaller c 10 1. A natural candidate is the unzipping oc- 
curring for the diblock copolymer model studied in Ref. 
p3| . The physics of that system can be assimilated to 
that of a DNA molecule in which each of the two strands 
is made exclusively of one type of base, and the bases 
of the two strands are complementary Of course, 

in such a model the loops forming can be made with 
portions of different lengths of the two strands. More- 
over, different loops can also bind and form more com- 
plicated topologies than in the DNA case. By apply- 
ing the simple homogeneous network picture to such a 
model, one would expect a second order transition with 
c = di/ - 2(73 - 1 « 1.4 in d = 2. The results for 
the copolymer model of Ref. , suggest a slightly, but 
definitely higher value c = 1 + 9/16 « 1.56, which is 
supported by a connection with percolation theory. This 
discrepancy could be due to the more complicated topol- 
ogy of the loops and could indicate that the model is less 
favorable for the application of a network picture. 

Homogeneous networks are very interesting per se and 
most recently were recognized as important tools also 
for the study of the topological entanglement of poly- 
mers ||2^ . Block copolymer networks are a still relatively 
unexplored subject, in spite of the obvious fundamental 
and applicative interest. Thus, the realization that such 
networks could also be very relevant for a description of 
system like DNA at denaturation, at the moment, can 
not yield quantitative predictions based on field theory 
results. Indeed, even the problem of identifying which 
kind of copolymer model in the continuum, if any, could 
represent correctly our discrete DNA in the scaling limit 
is far from trivial and completely open. For sure the 
identification of DNA as a system potentially connected 
to the copolymer network physics adds further interest 
to these intriguing objects, which are already known to 
be somehow related to multifractal aspects of polymer 
statistics and, in the d = 2 case, could realize inter- 
esting examples of conformal invariance. 

We thank D. Mukamel for fruitful discussions and 
E. Orlandini for ongoing collaboration. Financial sup- 
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PAIS 2001 and by European Network on "Fractal and 
Structures and Self-Organization" are gratefully ac- 
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